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^ : ABSTRACT 
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\^ ' We have analyzed the profile of the millisecond pulsar PSR J0437-4715 at 

I 1440 MHz by fitting the Gaussians to pulse components, and identified its 11 

^ ' emission components. We propose that they form a emission beam with 5 nested 

Q"^ ' cones centered on the core. Using the phase location of component peaks, we 

J^ ■ have estimated the aberration-retardation (A/R) phase shift. Due to A/R phase 

Tij- ■ shift, the centroid of intensity profile and the inflection point of polarization 

^ ' angle swing are symmetrically shifted in the opposite directions with respect to 

O ■ the meridional plane, which is defined by the rotation and magnetic axes. By 

1-^ ■ recognizing this fact, we have been able to locate the phase location of meridional 

I ' plane and estimate the absolute altitude of emission of core and conal components 

;_i ■ relative to the neutron star center. Using the more exact expression for phase 

^ ■ shift given recently by Gangadhara (2005), we find that the radio emission comes 

IJ • from a range of altitude starting from the core at 7% of light cylinder radius to 

outer most cone at 30%. 
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Subject headings: pulsars: general — pulsars: individual (PSR J0437-4715) 



1. INTRODUCTION 

The shape and polarization of integrated pulse of a pulsar are quite stable, and they 
are expected to the reflect the magnetic field geometry and the structure of pulsar emission 
beam. The common occurrence of odd number of components and their distribution across 
the pulse window has led to the nomenclature of core and cone structure for the pulsar 
emission beam (Rankin 1983a). However, Lyne & Manchester (1988) and Manchester (1995) 
have argued for a different interpretation, in which components on either side of core arise 
from a random distribution of emitting regions, which lead to the so-called patchy cone. 
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Mitra & Deshpande (1999) have attempted to test these models, and found the evidences 
in favor of conal structure. The frequency evolution of profile structure and polarization 
angle swing of pulsars has been interpreted by considering the radio emission by relativistic 
plasma moving in dipolar magnetic field lines (e.g., Radhakrishan & Cooke 1969; Sturrock 
1971; Ruderman & Sutherland 1975; Rankin 1983a&b, 1990, 1993; Lyne & Manchester 1988; 
Blaskiewicz et al. 1991; Hibschman & Arons 2001; Mitra & Li 2004). 

For understanding the pulsar radio emission mechanism, the radio emission altitude 
bears an at most importance. There have been two kinds of methods proposed for estimating 
the radio emission altitudes: (1) Purely geometric method, which assumes the pulse edge 
is emitted from the last open field lines (e.g., Cordes 1978; Gil & Kijak 1993; Kijak & 
Gil 2003), (2) Relativistic phase shift method, which assumes the asymmetry in the conal 
components phase location, relative to core, is due to the aberration-retardation phase shift 
(e.g., Gangadhara & Gupta 2001, hereafter GGOl). Both the methods have merits and 
demerits: former method has an ambiguity in identifying the last open field lines while the 
later rely on the clear identification of core-cone structure of the emission beam. Thomas & 
Gangadhara (2005) have showed that the profiles becomes asymmetric due to the infiuence 
of rotation on the radio emitting sources moving along the dipolar magnetic field lines. 

By assuming the pulsar radio emission comes from fixed altitude, Blaskiewicz et al. (1991) 
have showed that polarization angle infiection point (PAIP) lags with respect to the centroid 
of the intensity profile by ~ 4r/rLc, where r is the emission altitude measured from the 
center of neutron star and tlc = Pc/27r is the light cylinder radius. The parameters P and 
c are the pulsar spin period and speed of light, respectively. The results of purely geometric 
method are found to be in rough agreement with those of Blaskiewicz et al. (1991). How- 
ever, the relativistic phase shift method clearly indicates that the emission altitude across 
the pulse window is not constant (GGOl; Gupta & Gangadhara 2003, hereafter GG03; John- 
ston & Weisberg 2006). Dyks, Rudak and Harding (2004, hereafter DRH04) by revising the 
relation for aberration phase shift given in GGOl, have re-estimated the emission heights. 
The average emission heights are found to be comparable with those of geometric method. 

By considering the relativistically beamed radio emission in the direction of magnetic 
field line tangents, Gangadhara (2004, hereafter G04) has solved the viewing geometry in an 
inclined and rotating dipole magnetic field. He has showed that due to geometric restrictions 
in receiving the beamed emission, a distant observer can not receive the radio waves from 
a constant altitude. A more exact expression for relativistic phase shift (A/R) is given by 
Gangadhara (2005; hereafter G05), which includes the phase shift due to polar cap currents 
also. Under small angle approximations, the expression of G05 reduces to that of DRH04. 

Recently, Johnston & Weisberg (2006) have estimated the emission heights of the young 
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pulsar PSR J1015-5719 using the expression for A/R phase shift given in G05. In this 
paper, we intend to estimate the emission heights of the nearby, bright milhsecond pulsar 
PSR J0437-4715. It was discovered in the Parkes southern survey (Johnston et al. 1993; 
Manchester et al. 1996), has a pulse period of 5.75 ms. Navarro et al. (1997) have presented 
the mean pulse observations of PSR J0437-4715 at multiple radio frequencies, and found a 
frequency dependent lag of the total-intensity profile with respect to the polarization profile. 
Jenet et al. (1998) have presented the results based on the high time resolution single pulse 
observations of this pulsar. They find the single pulse properties of PSR J0437-4715 are 
similar to those of the common slow-rotating pulsars. 

The outline of this paper is as follows. In §2, we discuss a method for identifying 
the phase location of meridional plane. Next, in §3 we consider the mean pulse data of 
millisecond PSR J0437-4715, and fit the Gaussian to pulse components. By assuming conal 
emission beam, we have estimated the A/R phase shifts and the emission heights. 



2. Phase of meridional plane 

In the nested cone structure of pulsar emission beam, the conal components are expected 
to be distributed symmetrically on either side of core when observed in the corotating frame. 
However, for an observer in the laboratory frame the profile becomes asymmetric due to 
rotation. A typical asymmetric profile is given in Figure 1, and such an asymmetry in the 
real profiles has been reported for many pulsars (e.g., GGOl; GG03; Johnston & Weisberg 
2006). 

Blaskiewicz et al. (1991) have generalized the rotating vector model (RVM) by consid- 
ering the first order relativistic effects, which arise due to the pulsar rotation. By assuming 
that the emission across the entire pulse comes from a constant height, they have showed 
that the centroid of the total intensity pulse advances to earlier phase by ~ t/ti^q while the 
polarization angle inflection point (PAIP) is delayed to later phase by ~ 3r/rLc, where r 
is the emission height measured from the center of neutron star. However, there are clear 
evidences to show that the emission do not come from a fixed altitude (GGOl; GG03; John- 
ston & Weisberg 2006). Consider the intensity profile given in Figure 1, in which, Ci and 
Ct are the components of an arbitrary cone with center at CC. If the emission altitude is 
not constant across the pulse window then the retardation also contributes to the phase 
shift, and it is comparable to that of aberration (G05). So, due to A/R effects, the centroid 
intensity profile is advanced to earlier phase by ~ 2r/rLc while the corresponding PAIP 
is delayed to later phase by ~ 2r/rLc (e.g., DRH04; Fung 2005). Hence, the cone center 
CC and PAIP lie symmetrically on either side of meridional plane M, which is at the phase 



-4- 



0' = 0°. The meridional plane is defined by the plane containing the rotation axis and the 
magnetic axis. Note that the phase location of M is invariant with respect to r. In other 
words, both CC and PAIP come closer to M at smaller r, and move away from it at larger r. 
In the co-rotating frame, i.e., in the absence of rotation, the A/R effects vanish, and hence 
core peak and PAIP appear at the meridional plane M. 

The field aligned polar cap current do not introduce any significant phase shift into the 
phase of PAIP but it does introduce a positive offset into polarization angle, however, it 
roughly cancels due to the negative offset by aberration (Hibschman & Arons 2001). The 
phase shift of pulse components due to polar cap current is estimated recently by G05, and 
found to be quite small compared to A/R phase shift. 



3. Emission height in millisecond pulsar: PSR J0437-4715 

Consider the nearest bright millisecond pulsar PSR J0457-4715, which has a period of 
5.75 ms. It was discovered in Parkes southern survey by Johnston et al. (1993), and is 
in a binary system with white dwarf companion. Manchester & Johnston (1995) have pre- 
sented the mean pulse polarization properties of PSR J0437-4715 at 1440 MHz. There is 
significant linear and circular polarization across the pulse with rapid changes near the pulse 
center. The position angle has a complex swing across the pulse, which is not well fitted with 
the rotating vector model, probably, due to the presence of orthogonal polarization modes. 
The mean intensity pulse has a strong peak near the center, where the circular polarization 
shows a clear sense reversal and the polarization angle has a rapid sweep. These two features 
strongly indicate that it is a core. The profile shows more than 8 identifiable components. 



Table 1. Parameters related to radio emission from PSR J0437— 4715 at 1440 MHz 
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'Emission heights eoniputed from the approximate formula of DRH04, and ''those using the more exaet formula of G05. 
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Consider the average pulse profile given in Figure 2 for PSR J0437-4715 at 1440 MHz. 
To identify pulse components and to estimate their peak locations we followed the procedure 
of Gaussian fitting to pulse components. Two approaches have been developed and used 
by different authors. Unlike Kramer et al. (1994), who fitted the sum of Gaussians to the 
total pulse profile, we separately fitted a single Gaussian to each of the pulse components. 
We used the package Statistics' NonlinearFit' in Mathematica for fitting Gaussians to pulse 
components. It can fit the data to the model with the named variables and parameters, 
return the model evaluated at the parameter estimates achieving the leastsquares fit. The 
steps are as follows: (i) Fitted a Gaussian to the core component (VI), (ii) Subtracted the 
core fitted Gaussian from the data (raw), (iii) The residual data was then fitted for the next 
strongest peak, i.e., the component IV, (iv) Added the two Gaussians and subtracted from 
the raw data, (v) Next, fitted a Gaussian to the strongest peak (IX) in the residual data. 
The procedure was repeated for other peaks till the residual data has no prominent peak 
above the off pulse noise level. 

By this procedure we have been able to identify 11 of its emission components: I, II, 
III, IV, V, VI, VII, VIII, IX, X and XI, as indicated by the 11 Gaussians in Figure 2. The 
distribution of conal components about the core (VI) reflect the core-cone structure of the 
emission beam. We propose that they can be paired into 5 nested cones with core at the 
center. 

Manchester & Johnston (1995) have fltted the rotating vector model (RVM) of Rad- 
hakrishnan & Cooke (1969) to polarization angle data of PSR J0437-4715 and estimated the 
polarization parameters a = 145° and ( = 140°. The parameter a is the magnetic axis incli- 
nation angle relative to the pulsar spin axis, and ( = a + f3, where (3 is the sight line impact 
angle relative to the magnetic axis. However, for these values of a and ( the colatitude s/sid 
(see eq. [15] in GGOl) exceeds 1 for all the cones. The parameter s gives the foot location 
of field lines, which are associated with the conal emissions, relative to magnetic axis on the 
polar cap. It is normalized with the colatitude sm of last open field line. However, Gil and 
Krawczyk (1997) have reported a = 20° and [3 = —4° by fitting the average intensity pulse 
profile rather than by formally fitting the RVM. Further, they find the relativistic RVM by 
Blaskiewicz et al. (1991) calculated with these a and f3 seems to fit the observed position 
angle quite well than the non-relativistic RVM. Further, Gil and Krawczyk (1997) have re- 
ported that the centroid of intensity profile and PAIP are separated by ~ 40°. By taking the 
cone center (CC) of outermost cone as the centroid of intensity profile, we can locate the 
phase of meridional plane M. It is at the mid point, i.e., at 6(f)' ~ 20° between CC and PAIP. 
Using 6(f)' ~ 2r/ri^c, we get the emission height r ^ 48 km. We have adopted these values 
of a and (3 in our model as they also confirm s/s\of < 1 (see below). The phase 0' = 0°, 
in Figure 2, corresponds to the phase of meridional plane (M), which is at the mid point 
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between CC of outermost cone and PAIP of relativistic polarization angle swing. 

In column 1 of Table 1 we have given the cone numbers. The cone number stands for 
core component. The peak locations of conal components on leading and trailing sides are 
given in columns 2 and 3, respectively. The conal components are believed to arise from the 
nested hollow cones of emission (Rankin 1983a, b, 1990, 1993), which along with the central 
core emission, make up the pulsar emission beam. Let 0; and 0J be the peak locations of 
conal components on leading and trailing sides of a pulse profile, respectively. Then, using 
the following equations, we estimate the phase shift ((50') of cone center with respect to M 
and the phase location {</)') of component peaks in the absence of phase shift, i.e., in the 
corotating frame (see eq. [11] in GGOl): 

s<i>' = l{<P't + <p'i) , 0' = ^(0;-0;). (1) 

In column 4 of Table 1, we have given the values of 6(f)'. It increases in magnitude from inner 
most cone to 4th cone, but for outermost cone i.e., for 5th cone it decreases slightly. We 
suspect this decrease in 6(f)' for the 5th cone is due to the fact that emission spot moves 
closer to the rotation axis, and hence aberration phase shift becomes smaller (G05). 

The half-opening angle F (see eq. [7] in G04) of the emission beam is given in column 5, 
and it is about 30° for the last cone. In column 6, we have given the emission heights 
computed from the approximate expression given in DRH04. Using the more exact formalism 
for A/R phase shift given in G05 (see eq. [45]), we computed the emission heights and given 
in column 7. Their percentage values in tlc are given in column 8. It shows the emission 
in PSR J0437-4715 occurs over a range of altitude starting from the core at 20 Km to the 
outermost cone at 90 Km. In coulmn 9, we have given the refinement A(%) = (1 — 05/07)100, 
where cq and cj are the emission heights given in columns 6 and 7, respectively. It shows the 
refinement is quite significant for the outer cones. In column 10, we have given the colatitude 
s/siof of foot field lines, which are associated with the component emissions. It shows due to 
the relativistic beaming and geometric restrictions, observer tends to receive the emissions 
from the open field lines which are located in the colatitude range of 0.2 to 0.6 on the polar 
cap. 



4. DISCUSSION 

A long standing question in pulsar astronomy has been the location of the radio emission 
region in magnetosphere. The emission heights of PSR B0329+54 given in GGOl, six other 
pulsars in GG03 and the revised ones in DRII04 are all relative to the emission height of 
core. They assumed that the core emission height is negligible compared to the light cylinder 
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radius of classical pulsars, but it may not be so in the case of millisecond pulsars. Though 
Blaskiewicz et al. (1991) have assumed a constant emission height in their model, they have 
indeed indicated the possibility of core emission at a lower altitude than conal emissions. 
We identify, due to A/R phase shift, that the meridional plane (M) is located at the mid 
point between the phase of core peak and polarization angle inflection point (PAIP). By 
recognizing this fact, we have been able to estimate the absolute emission heights of core 
as well as conal components. We find that the core emission in PSR J0437— 4715 is indeed 
come from a lower altitude than the conal emissions. Note that our estimates of error bar 
do not include the error bars in a and f3, as they are not reported. For comparison we have 
reported the emission heights (see column 6 in Table 1) estimated from the approximate 
relation of DRH04, which is independent of the parameters a and j3. In the case of core and 
inner cones (1 and 2), the refinement in the emission height is < 12% (see column 9). 

We selected PSR J0437-4715 because of its strong core and clear conal emissions. Po- 
larization angle (PA) is complicated near the core component but we can still locate the 
meridional plane by avoiding fitting that part of PA data. The idea is that without the 
knowledge of polarization angle data near the core it is possible to locate the meridional 
plane by just fitting the relativistic polarization angle formula of BCW91 to segments of 
PA angle data corresponding to the conal components. Since the emission altitude for each 
cone is different, one can not fit the entire PA data with a constant altitude. If we can best 
fit a segment of PA data corresponding to a outer cone or any other inner cone it is possi- 
ble to locate the phase of the meridional plane. Since we used the formula of Blaskiewicz 
et al. (1991), in which the aberration phase shift is independent of phase, our results are 
approximate. 

The polarization angle curve obtained by fitting the segments of PA data under C/ and 
Ct gives a virtual infiection point (PAIP) corresponding certain height reported by fitting 
routine. Since the emission height vary for each cone, we shall have as many number of 
virtual infiection points as the cones. The point to be noted is that the cone center and 
PA virtual infiection point of each cone is approximately symmetrically located on either 
side of the meridional plane. Since the cone centers and the virtual PA infiection points 
are approximately symmetrically located on either side of meridional plane, it is possible to 
locate the phase of meridional plane even if we are unable to detect virtual PA infiection 
point for certain cones or core. 

The geometrical method, which by comparing the measured pulse widths with geomet- 
rical predictions from dipolar models, can yield absolute emission heights. But such absolute 
values are likely to be underestimated due to an ambiguity in identifying the last open field 
lines. The profile widths generally measured at a 10%-level can give the lower limits for 



emission altitudes (Xilouris et al. 1996). In PSR J0437— 4715 the components I and XI are 
at the intensity level of 1 to 2%, and the field lines associated with them are not last open 
field lines, but they are found to be located at the colatitude s/siof ~ 0.6 on polar cap. 

It is possible that the A/R phase shift can be reduced by the rotational distortion of 
magnetic field sweep back of the vacuum dipole magnetic field lines. It was first considered in 
detail by Shitov (1983). Dyks and Harding (2004) have investigated the rotational distortion 
of pulsar magnetic field by making the approximation of a vacuum magnetosphere. For 
0' = 110°, (3 = —4° and a = 20° we computed the phase shift due to magnetic field 
sweep back Scf)'^^^^^ (see eq. [49] in G05). But it is found to be < 0.6° for r/rLc < 0.3, 
which is much smaller than the aberration, retardation and polar cap current phase shifts 
in PSR J0437-4715. 



5. CONCLUSION 

We have analyzed the mean profile of the millisecond PSR J0437-4715 at 1440 MHz, 
and identified 11 of its emission components. We propose that they form a emission beam 
consist of 5 nested cones centered on the core. Using the phase location of component peaks, 
we have estimated the aberration-retardation (A/R) phase shift. Due to A/R phase shift, 
the center of each cone and the inflection point of corresponding polarization angle swing 
are symmetrically shifted in the opposite directions with respect to the meridional plane. By 
recognizing this fact, we have been able to locate the phase location of meridional plane and 
estimate the absolute altitudes of emission of core and conal components relative to neutron 
star center. We used the more exact expression for phase shift given recently by G05. The 
radio emission at 1440 MHz in PSR J0437-4715 is found to come from an altitude range 
starting from the core at 7% of light cylinder radius to the outer most cone at 30%. 

We used the data available on EPN archive maintained by MPIfR, Bonn. We are 
thankful to M. N. Manchester and S. Johnston for making their data available on web. 
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Fig. 1. — A typical pulse profile. Cq is the core. Ci and Ct are the components on leading 
and trailing sides, respectively. The dotted lines represent the phase of cone center (CC) 
and polarization angle inflection point (PAIP). The thick line mark the phase of meridional 
plane (M), and the broken line mark the phase of core peak (C). Arrows indicate direction 
of phase shift of CC, C and PAIP with respect to M. 
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Fig. 2.— Average pulse of PSR J0437-4715 at 1440 MHz and the model Gaussians (dotted 
line curves) fitted to emission components. Due to A/R effects, the core (VI) peak is shifted 
to the phase ~ —9.5° from the meridional plane (M) which is at 0°. 



